Kinetics of step bunching during growth: A minimal model 
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We study a minimal stochastic model of step bunching during growth on a one-dimensional 
vicinal surface. The formation of bunches is controlled by the preferential attachment of atoms to 
descending steps (inverse Ehrlich-Schwoebel effect) and the ratio d of the attachment rate to the 
terrace diffusion coefficient. For generic parameters (d > 0) the model exhibits a very slow crossover 
to a nontrivial asymptotic coarsening exponent f3 ~ 0.38. In the limit of infinitely fast terrace 
diffusion (d — 0) linear coarsening (/3 = 1) is observed instead. The different coarsening behaviors 
are related to the fact that bunches attain a finite speed in the limit of large size when d = 0, 
whereas the speed vanishes with increasing size when d > 0. For d = an analytic description of 
the speed and profile of stationary bunches is developed. 

PACS numbers: 81.10.-h, 05.40.-a, 68.55.-a, 89. 75. Da 



I. INTRODUCTION 

Step bunching is a morphological instability of a vic- 
inal crystal surface in which a regular train of equally 
spaced steps decomposes into alternating low step den- 
sity regions (terraces with orientation close to a singular 
surface) and high step density regions (bunches) [l| . Dur- 
ing the evolution the surface self-organizes into a pattern 
with a characteristic length scale, which coarsens in time. 
The phenomenon of considerable interest for applications 
|2| as well as from the fundamental point of view of non- 
equilibrium statistical mechanics. 

Step bunching is generally caused by breaking the sym- 
metry between the ascending (upper) and descending 
(lower) steps bordering a terrace. This can be due to 
a variety of mechanisms. In growth or sublimation, the 
symmetry is broken by different kinetic rates for the at- 
tachment or detachment of particles at the upper and the 
lower step Q, while in electromigration Q the asymme- 
try is introduced by the electric field. Other mechanisms 
responsible for step bunching during growth include im- 
purities 0, , diffusion anisotropy Q , and the presence 
of a second surface species @ . Quite generally, a growing 
equidistant step train shows a step bunching instability 
when adatoms attach to the step more easily from the 
upper terrace than from the lower terrace (an inverse 
Ehrlich-Schwoebel effect, iES). While iES behavior is dif- 
ficult to justify microscopically, it may serve as a use- 
ful effective description of more complex step bunching 
mechanisms 0. 

The coarsening of the bunched surface can be charac- 
terized by the power law increase of the average distance 
between bunches C as C oc t 13 . A few experimental ob- 
servations of step bunch coarsening during growth have 
been reported 0, 0, El 0] , but quantitative results are 
scarce. Numerical simulations of the equations of step 
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FIG. 1: Illustration of the model. Particles are deposited 
onto the terraces and attach to the ascending (descending) 
step with probability p+ (P-)- The attachment probabilities 
depend on the terrace width I through 



motion found (3=1/2 independent of the step interaction 
potential |l2j |. Recently a unifying continuum treatment 
of step bunching instabilities was proposed pH , FLU HE) 
with the aim to identify universality classes. Neverthe- 
less, a good understanding of the relationship between 
microscopic step bunching mechanisms and the resulting 
asymptotic scaling properties is so far lacking. 

In this paper, we study coarsening of step bunches dur- 
ing growth for a simple one-dimensional lattice model. 
The main idealization compared to conventional one- 
dimensional models [12( is that repulsive step-step inter- 
actions are ignored, and hence the individual steps can 
coalesce to form composite steps. This eliminates the 
internal structure of the bunch and the associated addi- 
tional length scale [l4|, and allows us to simulate large 
systems over very long times. Related models have been 
proposed previously [jj, fToL ll6L Il7| , but their asymptotic 
bunching behavior has not been explored. 



II. MODEL 

We consider S steps of unit height located at positions 
afi, X2, xs of a one-dimensional lattice with L sites 
and periodic boundary conditions; the slope of the vici- 
nal surface is mo = S/L. The width of the z-th terrace 
between steps i and i — 1 Due to the ab- 



2 



sence of step-step repulsion two or more steps can be lo- 
cated at the same position; correspondingly some terrace 
widths may be zero. Particles are uniformly deposited 
onto the sample. A particle deposited on a terrace of 
width I moves immediately to one of the bordering steps 
and is incorporated with the probabilities 



P± 



ll±b+d 
2 l + d£ 



(1) 



at the ascending (p+) or descending (p_) step, respec- 
tively (FigH)l. The balance b is a measure for the attach- 
ment asymmetry, while the inverse diffusivity d deter- 
mines how the asymmetry depends on the width of the 
terrace: For d = the attachment probabilities p± are 
independent of the terrace width, while for d > the at- 
tachment becomes symmetric when £ ^> 1/d. Note that 
in the symmetric case (b = 0) the value of d is irrelevant. 
For d = and 6=1 the model is exactly solvable, and 
the terrace widths have a Poisson distribution [l^ ■ 

The form (JIJ follows from the solution of the statio nary 
diffusion equation for the adatoms on the terrace 0, Il8j |. 
which leads to the following expressions for the model 
parameters in terms of the adatom diffusion coefficient 
D and the attachment rates k± 0,0 



k + — k- 



d : 



D/k + + D/k^ 



(2) 



(3) 



A normal Schwoebel effect (k + > fc_) corresponds to 
b > 0, while an iES effect implies b < 0. The quantity d 
is the inverse of the sum of the kinetic lengths D/k± [l| . 

As initial condition, the steps are placed at random. 
We use sample sizes ranging from L = 5 x 10 4 to 10 6 . 
Time is measured in number of monolayers (ML), and 
simulations were performed up to 10 7 ML. As expected, 
the model exhibits step bunching when b < 0. An exam- 
ple of the formation and time evolution of bunches for 
d = is shown in Fig. [21 Bunches and steps move to the 
right. Eventually in the stationary regime only a single 
bunch remains which moves with a constant velocity. In 
the case d > the coarsening is slower but the space- 
time plot is qualitatively similar. We observed that the 
bunch velocity is generally increasing with decreasing d 
and increasing \b\ (see Sect II VI for further discussion). 

The simplest way to characterize the surface morphol- 
ogy is to measure the surface width W relative to the 
mean tilted surface 0] . For 6 > no step bunches form 
and the surface is roughened only by fluctuations. For 
b > the model belongs to the Edwards- Wilkinson uni- 
versality class and W ~ t 1 / 4 ,1?]. For b — we find nu- 
merically that W ~ t 1 / 3 , which is consistent with the idea 
[20j that the symmetric model belongs to the conserved 
Kardar-Parisi-Zhang universality class 0|. In the step 
bunching regime b < the surface width should asymp- 
totically become proportional to the bunch spacing C. 
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FIG. 2: Spatiotemporal evolution of the step configuration on 
a vicinal surface of mean slope mo = 0.2. Time is increasing 
from bottom to top. The positions of the steps are shown as 
dots. The sample size in the horizontal direction is L = 5000. 
Parameters are d = and b = —0.6. 
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FIG. 3: Squared surface width W 2 and maximal terrace size 
Cax as a function of time for a system of size L — 2 x 10 5 , 
slope 0.1, and parameters b — — 1, d = 0.1. The bold lines 
have slope 2/3 = 0.76 and (3/2 = 0.19, respectively. 



Indeed we find that W increases faster than t 1 / 3 when 
b < (see Fig. |3J), but the large background contribu- 
tion to W which arises from random fluctuations makes 
this quantity less suited for a precise determination of 
the coarsening behavior. We therefore developed an al- 
ternative coarsening measure which we describe next. 
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FIG. 4: Density-density correlation function for a system of 
size L = 3000, parameters b = —0.3 and d = 0.01, and mean 
slope of the vicinal surface equal to 0.1. The data are averaged 
over 500 independent runs. 



III. COARSENING BEHAVIOR 

There arc different possible ways to quantitatively ana- 
lyze the coarsening process. In analogy with phase order- 
ing kinetics [2lJ , we calculated the density-density corre- 
lation function 



C n (x,t) = j y^(n(y, t)n(y + x, t)) , 



L 

v 

the step density being defined as 

S 



(4) 



i(x,t) = ) 6(x- Xi{t)) 



(5) 



i=l 



where 8(x) — 1 for x — and zero otherwise. 

The typical time evolution of the correlation function is 
plotted in Fig. We observed a minimum in the density- 
density correlation function, which shifts to larger values 



during growth and ultimately saturates at half the sys- 
tem size. We did not observe any significant secondary 
maximum in C n {x 1 1), indicating that there is no system- 
atic periodicity in the bunch configurations. The position 
of the minimum is a good indicator of the typical bunch 
distance. On the other hand, it is hardly measurable 
with sufficient precision. Even after averaging over 500 
independent runs the results are not satisfactory. There- 
fore, we turned to other quantities to effectively capture 
the dynamics of bunching. 

We may operationally define a composite step of size 
k as an object in which there are at least k steps at the 
same position. For k = 1 we are dealing with the steps 
themselves, while k > 1 corresponds to composite steps 
in the strict sense. Denote by Sk the number of such 
composite steps and Xk,i,i — 1,2, ...,Sk their positions. 
If we identified the step bunches with composite steps 
and measured the average bunch distance as L/Sk, we 
would make a systematic error, because often two such 
"bunches" come very close one to the other, effectively 
making up a single larger bunch. We need to eliminate 
these short distances. To this end we define the quantity 



A fe = 2^,(( x k,i 
i=l 



Xk,: 



-if), 



(6) 



in which larger distances dominate over shorter ones. 
Therefore, the quantity Ai measures the effective dis- 
tance between steps in regions of low density of steps, 
i.e. around the middle between the bunches. On the 
other hand, A& for k > 1 measures the effective distance 
between bunches, irrespective of the elementary steps dis- 
persed between bunches. 

We compared the time evolution of for various val- 
ues of k. As expected, the behavior differs qualitatively 
for k = 1 and for k > 1. The quantity Ai grows with 
time and eventually saturates at a value which is signif- 
icantly lower than the system size. The saturated value 
increases with increasing \b\ and approaches the system 
size for b — * — 1. 

On the other hand, the quantities A& for k > 1 exhibit 
first a transient period with a fast decrease and then grow 
until they reach a value very close to the system size. (In 
principle the quantity is always lower than the system 
size, but we observed that this difference becomes negli- 
gible when we increase the system size. This fact can be 
easily undersood from the form of the stationary bunch 
profile, as will be shown below in Sect II V B)| . The be- 
havior depends on k only in the initial transient period, 
which is longer for larger k. However, the universal long 
time behavior, which is the subject of this work, is found 
to be independent of k > 1. Therefore, we will concen- 
trate on the quantity A2 in the following. 

We investigated the time evolution of the distance be- 
tween bunches in different regimes depending on the pa- 
rameters b and d. Typical results are shown in Fig. 
for b = —0.9 and several values of the diffusion parame- 
ter d (the results for other values of b < are similar). 
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FIG. 5: Time evolution of the distance between bunches mea- 
sured by A2 for b — —0.9, and varying diffusion parameter 
d = 0, 0.0001, 0.001, 0.003, 0.01, 0.1, and 1 (from top to bot- 
tom at time 10 5 ). System sizes are 10 6 , 10 6 , 5 x 10 5 , 5 x 10 5 , 
10 5 , 10 5 , and 5 x 10 4 , respectively. The lines are powers cx t 1 
and t ' 38 . The mean slope of the vicinal surface is 0.2. 

For d — 0, linear behavior A2 oc t is reached after some 
transient, i.e. (3 ~ 1. On the other hand, the situa- 
tion for nonzero d is more complex. After the initial 
transient, there is an intermediate behavior with approx- 
imately power-law character, but the exponent strongly 
depends on the value of d, while the dependence on b was 
very weak. 

Eventually, after some crossover time ~ t c a stationary 
power-law regime is reached, which in our data extends 
over several orders of magnitude in time. The effective 
exponent f3 is only weakly varying with the nonzero d and 
assumes the asymptotic value (3 ~ 0.38. We conjecture 
that this regime is universal and the remaining variations 
in the value of the exponent are solely due to finite size 
effects. 

In order to confirm that the asymptotic behavior is uni- 
versal, one can rescale the raw data for different nonzero 
rf, plotting them as functions of t/t c in the form 

A 2 (t; b, d) = A 2z + (t/t c )°- 3S A 2c F(t/t c ). (7) 

The results for three values of d and three values of b 
are shown in Fig. Neglecting the initial transient, 
we see that the scaling funtion F(x) approaches unity 
for x ^> 1, marking the asymptotic universal regime, and 
behaves close to a power-law F(x) ~ x T for x <C 1, which 
characterizes the intermediate regime; the exponent r 
depends on d but is independent of b. The initial offset 
A 2z was introduced for convenience. 

The crossover time t c between the intermediate and 
asymptotic regimes turned out to be rather large (t c ~ 
10 4 for d = 0.001 and b = —0.9) and it increases rapidly 
with decreasing \b\, as can be seen in the inset of Fig. [U 
where behavior consistent with a power law t c ~ (— 6)~ 4 



FIG. 6: Scaled time evolution of the effective distance between 
bunches measured by A2 for three values of the balance and 
three values of the inverse diffusivity. Groups of curves cor- 
respond to diffusion parameters d — 0.01, 0.001, 0.0001 from 
top to bottom. Each group consists of data for b — —0.9 (□, 
V, ▼), -0.6 (x, A, A), and -0.3 (+, 0, • ). The slope of the 
vicinal surface and the sample sizes are the same as in Fig. 
[3] The inset shows the values of the crossover time t c used in 
the main plot. The values of d are 10 -4 (+), 10 -3 (x), and 
0.01 (□). The line is a power law oc (—ft) -4 . 



is apparent. The dependence on d turns out to be non- 
monotonic: The crossover time decreases up to about d = 
0.01 and then increases again. This shows that different 
crossover mechanisms are acting for small and large d, 
respectively. 

Indeed, for very small d the dependence of the attach- 
ment asymmetry on the terrace size is felt only once the 
terraces become sufficiently large, and the system be- 
haves initially as if d = and (3 = 1; correspondingly 
the coarsening exponent approaches its asymptotic value 
from above in this case. Conversely, for large d the asym- 
metry in the attachment probabilities is always very 
small; consequently it takes a long time for the bunches 
to form, and the asymptotic coarsening exponent is ap- 
proached from below. 

The asymptotically linear behavior of the bunch size 
for d = 0, (3 rs 1, is well confirmed by measurements of 
the surface width W, and the data for d = 0.1 shown in 
Fig. [3]are consistent with the estimate (3 « 0.38 obtained 
from A2. In addition, Fig|3] displays the evolution of the 
maximal terrace width ^ max m the system. This length 
scale is seen to scale with a distinct exponent which is 
approximately given by 0/2 ~ 0.19. We will return to 
the behavior of l max below in Sect llVBl 



■5 



50400 

50350 - 

50300 -_ 

50250 - 

50200 _ 

50150 - 

50100 - 

50050 - 

50000 — 




V 



00 400 600 800 1000 1200 1400 1600 

X 

FIG. 7: Height profiles of a stationary bunch for d = 0, b — 
—0.1. The figure shows a surface with S = 160, L = 1600 
after deposition of 50000, 50100 and 50200 ML, respectively. 



IV. STATIONARY BUNCHES 

The final state of the system for any b < is found 
to be a single bunch propagating at constant velocity 
(FigCJ. The analysis of these stationary configurations 
provides important clues to the different coarsening be- 
haviors observed for d — and d > 0, respectively. 



A. Bunch speed 

Figure [H] shows numerical results for the stationary 
bunch speed v. For d — 0, v approaches a finite limiting 
value Voo with increasing system size; this is also evident 
from the space-time plot in Fig|21 The bunch speed is 
easily determined for the limiting case d = 0, b = — 1. fn 
this limit single steps cannot detach from the bunches, 
and the formation of composite steps is irreversible. The 
motion of a bunch is solely due to the atoms deposited 
onto the trailing terrace behind the bunch, which has 
length L in the stationary state. To move the bunch 
containing all S steps laterally by one lattice spacing, S 
atoms have to be deposited on this terrace; the bunch 
speed is therefore Voo — L/S — 1/mo- The limiting 
bunch speed decreases with decreasing |6| and vanishes 
as Voo ~ b 2 for 6^0 (see inset of FigJSJ. An ana- 
lytic explanation of this behavior will be given below in 
Sect llVTI 

For d > the bunch speed decreases indefinitely with 
increasing L. The simulation results shown in FigJ^jare 
consistent with a behavior 



v{L) ~ L~ v 



(8) 



with v w 1, but owing to the significant curvature of 
the data, the asymptotic value of v cannot be accurately 
determined. 

Following a simple scaling argument due to Chernov 
|22^ . we can try to relate the exponent v in (JSJ to the 
coarsening exponent (3. At a time when the typical bunch 
spacing is C, the bunch speeds are of order L~ v . Assum- 
ing that there is only a single scale in the problem, the 
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FIG. 8: Propagation speed v of stationary bunches. The main 
plot shows the dependence of v on the system size L, contrast- 
ing the behavior for d = and d > 0. The inset shows the 
bunch speed for d = 0, extrapolated to L = oo, as a function 
of \b\; the extrapolation is based on data up to L = 1600 and 
assumes a leading correction of order 1/L. All data were ob- 
tained by time averages in the stationary regime. The surface 
slope was 0.1 for d = and 0.2 for d = 0.01. The full bold 
line illustrates the behavior v ~ 1/L. The dashed bold lines 
are the upper bounds on v derived in Sect II V CI The bounds 
in the main figure were derived from the algebraic equation 
12 It , while the bound in the inset is the inequality I16H which 
becomes exact for 161 — > 0. 



velocity difference Aw between two bunches is of the same 
order. The time t* until the coalescence of two bunches 
can then be estimated as t* ~ £/ Av ~ and revers- 

ing this relationship one obtains the expression 



(9) 



for the coarsening exponent. For d = we have seen that 
the bunch speed remains finite for L — > oo, hence v = 
and P = 1 in agreement with the coarsening simulations. 

For d > the value v 1 implies (3 w 1/2, which is 
larger than our numerical estimate w 0.38. Although 
the data for v(L) in Fig|S]are probably not asymptotic, 
it seems unlikely that they will ever reach the large slope 
v « 1.6 required to reproduce our value for /3. In fact, 
our simulations show clearly that, besides the coalescence 
of bunches, the exchange of steps between bunches of 
different size plays an important role in the coarsening 
process. In addition, the existence of the distinct length 
scale < max (see FigOJ invalidates the assumption that the 
bunch spacing is the only length scale in the problem. 
Thus, while Chernov's relation (0 helps to connect the 
different coarsening behaviors for d = and d > with 
the difference in the size dependence of the bunch speed 
on a qualitative level, it is not quantitatively satisfied for 
our model. 
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FIG. 10: Maximal terrace size £ m ax for stationary bunches as 
a function of system size L. The mean slope is mo = 0.2. The 
data were obtained as time averages over single runs. 



FIG. 9: Stationary bunch profile at time 1 x 10 7 for sample 
size L = 2 x 10 J with slope 0.2 for d — and balance b — —0.3. 
Inset shows log-linear plot of profile relative to the center of 
bunch xo, averaged over 11 (for b = —0.3) and 5 (for b — —0.6) 
independent runs. The lines are predictions from the formulae 
itTUt and JSTt . 



B. Height profile 

We measured the height profile between bunches in the 
stationary regime. We show in Fig. [5]thc results for d = 
and two values b — —0.3 and b = —0.6. Qualitatively we 
observe that for larger values of —b the bunch becomes 
steeper in its very center, while the rest of the bunch, 
farther from the center, is more gradual. Since there 
is no repulsion between steps they tend to accumulate 
in the center of bunch (position denoted xq) and to the 
right of it, while the region on the left-hand side is much 
more flat. This means that the bunch shape is strongly 
asymmetric. Due to the discreteness of the lattice it it 
very difficult to infer precisely the singular behavior near 
the center of bunch, but the averaged data in the inset 
show that the profile far from the center can be described 
by a logarithmic behavior of the form 



h(xo) — h(x) » Ahi(x — xq). 



(10) 



The logarithmic singularity observed in the bunch shape 
is very weak and in all cases investigated most of the 
mass of the bunch is concentrated on two sites in its 
center. Indeed, setting x — Xq — L in l|10|) it is seen that 
the number St of terrace steps which are not contained in 
the center of the bunch grows only logarithmically with 
L. This is the ultimate reason why the effective quantity 
A2 measures quite well the distance between bunches. 

The logarithmic height profile implies that the surface 
slope m(x) — dh/dx decays as mix) rs —A/x far away 
from the bunch. The size f max of the largest terrace in 
the system, which is found at a distance of order L from 
the bunch, is then of order ^ max ~ —m(L) ~ L. This 
is consistent with measurements of Z max for stationary 



bunches with d = (Fig llOfl . On the other hand, for 
d > the largest terrace is much smaller, and scales 
approximately as € max ~ L 1 / 2 . This matches the time- 
dependent behavior of € max shown in FigOU 



C. Analytic description of stationary bunches 

In this section we provide an analytic derivation of sev- 
eral of the numerically observed properties of stationary 
bunches for the case d = [25l |. A convenient starting 
point is the continuum evolution equation 



dh d 
dt dx 



\b\ 1 dm 
2m 6m 3 dx 



(11) 



which has been obtained for the deterministic version of 
the present model throu gh a n essentially rigorous coarse 
graining procedure [Til l20|. Here m = dh/dx is the 
surface slope, the time scale is chosen such that the de- 
position flux equals unity, and b is assumed negative. A 
continuum description of the type (|1 1|) is expected to be 
valid far ahead of the bunch, where the step spacing is 
large compared to the atomic scale and lattice effects are 
negligible. 

We are looking for solutions of (|llf) which describe a 
bunch of height S (in units of the vertical lattice spacing) 
in a system of length L, moving laterally at speed v. This 
implies the ansatz 



h(x,t) = f{x - vt) + Qt, 



(12) 



where the last term accounts for the fact that also the re- 
gions between bunches grow vertically due to the terrace 
steps that move across these regions [23| (see FigJTJ . The 
constant unit flux on the right hand side of (|llf) implies 
the sum rule Hjl 



o S 

n + -v 

Li 



1 



(13) 



connecting the vertical and lateral growth rates. Insert- 
ing (|T2*|l into (fTTf> one obtains the ordinary differential 
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equation 



~ V 9 + 17 



\b[ 1 dg 
2g 6 5 3 dt; 



= i -n 



(14) 



where t; = x — vt and g = df /d£ is the surface slope in the 
comoving frame. The first term on the left hand side of 
(|14fl can be neglected relative to the right hand side when 
g <C mo, which is true far ahead of the bunch. Then it 
is readily verified that the equation admits a solution of 
the form g(£) = —A/^, consistent with the height profile 
QKljl. The coefficient A satisfies the quadratic equation 



{l-n)A 2 - A\b\/2 + 1/6 = 0. 



(15) 



Using JT3J, one finds that l|15f) has real solutions only 
when the bunch speed satisfies the inequality 



3 b 2 
v < . 

8 mo 



(16) 



Whereas the proportionality v ~ b 2 suggested by this 
relation is well confirmed by the data in the inset of FigJSJ 
the bound on the numerical coefficient seems to be weakly 
violated by the data. However, since the convergence of 
v with increasing L is quite slow, we believe that our 
simulations can still be regarded to be consistent with 
(|16fl . Assuming that l|16|) is satisfied essentially as an 
equality, we find that the coefficient A is given by 



2 



(17) 



The relations (|16|l and l|17fl can be improved by ana- 
lyzing the discrete evolution equations for the mean step 
positions X; = (xj). For d — they take the simple 
linear form |3ll20l| 

^ = i(l + &)(AVi-^) + ^(l-&)(^-^-i)- (18) 
The moving bunch is described by a solution of the form 



Xi(t) = <t>(i + Qt) + vt, 



(19) 



where f2 and v have the same meaning as in (H2J). Insert- 
ing (jT§|l into i(T%|) one obtains the difference-differential 
equation 

V + Q ^ = + + mC + 1} ~ \ {l bm ~ 1] 

(20) 

for the function <^(C)j where £ = i + fit. The expected 
logarithmic height profile (|10|1 corresponds to an expo- 
nential increase of <f> as (/)(() ~ e 1 */" 4 . For large Q the 
constant v on the left hand side of (|20|l can be neglected. 
Inserting the exponential ansatz into (121 )[) yields the al- 
gebraic equation 

n = A{smh(l/A) + |6|[1 - cosh(l/A)]}, (21) 



which reduces to (|15|l for small |6| and large A. For gen- 
eral values of |6|, l|21l) has real solutions only when the 
bunch speed v is smaller than an upper bound w m ax(&), 
which reduces to ()16J) for small \b\. In Figs |H1 and El we 
compare the predictions for v max and A derived from l|21|) 
with our numerical results. For A we use the value ob- 
tained from H21(l under the assumption that v = w max , 
which leads to a rather satisfactory agreement. 



V. CONCLUSIONS 

In this paper we have developed and studied a one- 
dimensional stochastic growth model which contains 
some essential features of the formation and coarsening 
of step bunches. Despite its simplicity, the model dis- 
plays two distinct scaling regimes. In the limit of fast 
terrace diffusion (d = 0) we have numerically established 
a linear coarsening law (/? = 1) and related this behavior 
to the fact that the propagation speed of bunches tends 
to a finite value with increasing bunch size. 

When the diffusivity is finite (d > 0) the attachment 
asymmetry decreases with increasing terrace size and the 
coarsening is slowed down. After having accounted for a 
complicated crossover behavior, we identified the asymp- 
totic value of the coarsening exponent, which we estimate 
to be (3 ~ 0.38. This exponent is different from exponents 
found in previous studies of step bunching, and it does 
not seem to follow from any simple scaling argument. It 
is robust, in the sense of being independent of the pre- 
cise values of the model parameters, provided b < and 
d > 0. 

Similar to earlier studies of kinetic roughening [26 ] , our 
work shows that finding the true asymptotic exponent 
may require a considerable amount of computer power 
due to a long intermediate regime with non-universal, ap- 
parent power-law behavior. The situation becomes even 
more difficult when we try to observe the crossover from 
j3 = 1 to f3 ~ 0.38, which should occur for sufficiently 
small but non-zero d and large enough system size L. 
Currently this is beyond our capacity. 

Due to the absence of step-step interactions, our model 
cannot provide an accurate description of the structure 
and dynamics of real step bunches. Nevertheless our 
study clarifies some conceptual issues that arise also for 
more realistic models. For example, although the inter- 
nal width of the bunches has been eliminated by allow- 
ing steps to coalesce, the model still displays (for d > 0) 
a second length scale, the maximum terrace size £ max , 
which grows with an exponent that is distinct from the 
coarsening exponent. This invalidates simple scaling ar- 
guments iKill 4l22i which assume the bunch spacing to 
be the only scale in the problem, and which would imply 
a " superuniveral" coarsening exponent (3 = 1/2 |27j. 

A second important advance of our work is the demon- 
stration that continuum evolution equations can be used 
to quantitatively describe the profile and speed of mov- 
ing step bunches. So far the use of continuum approaches 



8 



in step bunching has been limited to scaling arguments 
13] and to the study of stationary bunch configurations 

To conclude, we speculate that the relationship be- 
tween growing stepped surfaces and interacting particle 
systems (specifically, zero r ang e processes [2^) that was 
pointed out some time ago [2£| may be usefully exploited 
to gain further insight into the coarsening behavior of 
step bunches. In the language of interacting particle sys- 
tems, step bunching is an example of a condensation tran- 
sition, where a macroscopic fraction of the particles (i.e., 
the single height steps in our model) condenses into a 
few or a single site. The dynamics of such condensation 
transitions is of great current interest jlSE^I; and closely 



related concepts appear e.g . in granular physics [3(], EI 
and traffic research |2fJ, |32( . 
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